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Qh. Abstract 

Jr~^' We provide a comparative study of the fine tuning amount (A) at the two-loop leading 

r^ ■ log level in supersymmetric models commonly used in SUSY searches at the LHC. These 

are the constrained MSSM (CMSSM), non- universal Higgs masses models (NUHMl, 

psj I NUHM2), non- universal gaugino masses model (NUGM) and GUT related gaugino masses 

K*" ' models (NUGMd). Two definitions of the fine tuning are used, the first (Amax) measures 

vQ . maximal fine-tuning wrt individual parameters while the second (Ag) adds their con- 

l/^ I tribution in "quadrature" . As a direct consequence of two theoretical constraints (the 

EW minimum conditions), fine tuning (Ag) emerges at the mathematical level as a sup- 



CO . pressing factor (effective prior) of the averaged likelihood (£) under the priors, under 

^"^ the integral of the global probability of measuring the data (Baycsian evidence p{D)). 



For each model, there is little difference between Ag, Amax in the region allowed by the 
data, with similar behaviour as functions of the Higgs, gluino, stop mass or SUSY scale 
{fTisusy = ("^ti'TT'ts)^ ) '-'^ dark matter and g — 2 constraints. The analysis has the ad- 
vantage that by replacing any of these mass scales or constraints by their latest bounds 
^ _• one easily infers for each model the value of Ag, Amax or vice versa. For all models, 

minimal fine tuning is achieved for Mhiggs near 115 GeV with a Ag « Amax ~ 10 to 100 
depending on the model, and in the CMSSM this is actually a global minimum. Due to 
a strong (w exponential) dependence of A on Mhiggs, for a Higgs mass near 125 GeV, 
the above values of Ag « A,nax increase to between 500 and 1000. Possible corrections 
to these values are briefiy discussed. 



^E-mail: dumitru.ghilencea@cern.ch, hyun.min.lee@cern.ch, myeonghun.park@cern.ch 



1 Introduction 

Low energy (TeV-scale) supersymmetry (SUSY) can provide a solution to the hierarchy 
problem. This is done without undue amount of electroweak scale fine tuning (A) that 
is present in the non-super symmetric theories like the Standard Model (SM). A large value 
of this A is just another face of the hierarchy problem (for a review see [1] and references 
therein). However, negative searches for superpartners increase the SUSY scale (msusy) which 
in turn can increase A. In the extreme case when rususy is very high (>> TeV) one recovers 
the scenario of non-supersymmetric theories (SM, etc) with a large A. In the light of current 
negative SUSY searches at the LHC it is useful to examine in detail the amount of fine tuning 
that supersymmetric models need, as a test of SUSY as a solution to the hierarchy problem. 
The alternative is to ignore this problem and adopt an effective theory approach, with a low 
effective cutoff (few TeV) that, unlike SUSY, does not detail the "new physics" at/beyond 
this scale. Such models usually have a A relative to the TeV scale comparable to that of 
SUSY models relative to the Planck scale. 

While a small value of A (say less than 100) is desirable, the exact value still accepted for 
a solution to the hierarchy problem is rather subjective. Even worse, there are also different 
definitions of A in the literature. Two common definitions are 

Amaa; = max|A^J, A, = (^ ^ A^ J , A^^ = , 7^ = mo,mi/2,/io,^o, ^o- (1) 

7j are new parameters (of mass dimension 1), that SUSY introduces in the model (shown 
above for the CMSSM). A^ax was the first measure used [2], but Ag is also common. Two 
definitions for A can lead to different predictions and the absence of a widely accepted upper 
value for it is another problem. To avoid these issues, we compute both Ag and ^max and 
compare their implications in generic SUSY models, without an upper bound on them (to be 
fixed by the reader). This is one of the main purposes of this work. 

These measures of fine tuning were introduced more on physical intuition than rigorous 
mathematical grounds so another important purpose is to clarify their link with other ap- 
proaches and find technical support. Both A provide a local measure (in the space 7^) of the 
quantum cancellations and help to decide which phase space points of a model are less fine 
tuned (more probable). When actually comparing models, a more global measure would be 
desirable. Our scan over the whole parameter space when evaluating A's will alleviate this 
issue. But one question remains: what is the relation of A to other (global) measures of the 
success of SUSY in solving the hierarchy problem? To answer this, consider the Bayesian 
probability density V{'^i\D) of a point in parameter space {71} given the data D: 

Vijm = -^^C{D\^i) p(70, piD)= JciD\j,)pij,)dj,. (2) 

Here C is the likelihood the parameters {ji} fit the data D and p{D) is a global normalization 
factor called Bayesian "evidence". For two models with the same data and priors piji), the 
ratio of their p{D) gives their relative overall probability. So a large p{D) is needed to 



decide that a model is more probable than another. Then what is the relation between p{D) 
and fine-tuning? As it was observed in [^ (see also Uj), when integrating V{'ji\D) over one 
parameter of the theory (in this case fio), following an experimenta/ constraint (on mz), there 
is an emergent effective prior Pe// ~ l/'^fto which brings in a fine tuning penalty for points 
with largcl A^^ ~ A^ax 13]. These points then have little contribution to p{D) because 
J V ^ C X peff ~ (1/Amax) >C. A larger p{D) can then indicate a preference for points of 
lower Afnax, and the link oi p{D) with fine tuning wrt individual parameters is apparent. 

We explore this idea further and evaluate p{D) by investigating the effect of the theoretical 
constraints that received less attention: we refer to the two minimum conditions of the 
potential. Further, the above observation and the need to evaluate p{D) clearly suggests 
to integrate over all {7j} parameters (and we shall do so), and also over nuisance variables, 
which are parameters already present in the SM (like Yukawa couplings [31 [6]). The result is 
that p{D) ~ 1/Aq so Aq is actually preferred by the Bayesian evidence calculation. Then p{D) 
receives contributions mostly from points of small Ag, but this also depends on the priors 
and C To conclude, the inverse of Ag acts as an extra, effective prior in ([2]) and is indeed a 
physical quantity with impact on global p{D). This clarifies the exact, mathematical link of 
overall fine tuning wrt all parameters {7i}, to C{D\^i) and the Bayesian evidence p{D). 

With this technical motivation, we then evaluate the fine tuning for generic models, using 
both definitions Ag and A^ax] this is done in a two- loop leading log numerical analysis that 
provides the state-of-the-art analysis of the fine tuning in the models considered, consistent 
with current data. Interestingly, the results we find are little dependent on the definition used 
for A, with Ag and Amax showing similar behaviour and values. This is important since it is 
usually thought that different fine tuning measures should give different results. Our results 
correspond to a scan over the entire parameter space of the models (including tan/3). This 
is an extremely CPU-intensive task, made possible by the CERN batch computing service. 
The analysis of A^ax, Aq is done for the following models: 

a) - CMSSM: the constrained minimal supersymmetric standard model. For a recent 
two- loop leading log analysis of this topic see [5] , and for earlier investigations see [7] . 

b) - NUHMl: a CMSSM-like model but with equal Higgs soft masses, different from vriQ. 

c) - NUHM2: as for CMSSM but with different Higgs soft masses and different from niQ. 

d) - NUGM: a CMSSM-like model but with non-universal gaugino masses. 

e) - NUGMd: a benchmark NUGM model [8] with a GUT relation among gaugino masses. 

For these models our results are presented in a comparative way with Ag, Amax as func- 
tions of the lightest higgs, gluino, stop mass or the SUSY scale. Any experimental constraints 
on these can easily be used to identify Ag, Amax for that model. On top of these plots various 
contour lines corresponding to the remaining masses, dark matter or the g — l constraints are 
shown. Such comparative analysis for different models and definitions of A was not done in 
the past and has the advantage that it can be updated by the latest data, without re-doing 
the whole analysis. In particular, for each model we identify the corresponding A's for a 
Higgs mass of M^iggg = 125 ± 2 GeV that seems favoured by Atlas and CMS [U [10]. We shall 



^Note however that fine tuning wrt /lo is not dominant in CMSSM for higgs mass above '^ 115 GeV [5]- 



see such value requires Ag ?» ^max ~ 500 to 1000 depending on the model, and uncertainties 
in A are also discussed. In all cases A is minimal near M^igggK, 115 GeV. For ways to have 
Mhiggs-l'i^^ - 130 GeV with smaller fine-tuning A«C'(10) in SUSY models see [UlII]. 

In the following Section [2] shows the link of A to the evidence p{D) in models with theo- 
retical constraints. Numerical results and corresponding plots of A are shown in Section [3l 

2 Fine tuning, p{D) and the role of theoretical constraints. 

Before our numerical analysis, we re-examine the relation between the Bayesian probability 
of a point in parameter space or the evidence p{D) and the EW scale fine tuning, in models 
with theoretical constraints. Without loss of generality, we do this for the constrained MSSM 
(CMSSM). This analysis extends a previous similar study of this problem of [3] (Section 2 
in both papers), see also [3]. To place this discussion on quantitative grounds consider the 
CMSSM scalar potential 

,2 I rj |2 I ,y~^'2. I Tj |2 { ry-^'^ 



V = m{ \Hiy + m\ l^ar - im% Hi-H2 + h.c.) 



+ (Ai/2) \Hi\'^ + (A2/2) 1^2!^ + A3 \Hi\'^ \H2\'^ + A4 \Hi ■ H2'^ 



+ 



(A5/2) (Hi ■ H2Y + Ag \Hi\' {Hi ■ H2) + h \H2Y (Hi ■ H2) + h.c. . (3) 



The couplings \j and the soft masses receive one- and two-loop corrections that for the MSSM 
can be found in |121I13|. Let us introduce the notation 



m = rrii cos /3 -|- 777-2 sin /3 — m^ sin 2/3 
A = ^ cos^ /3 + ^ sin^ /? + ^ sin^ 2/3 + sin 2/3 {\ cos^ /3 + Ay sin^ /3) (4) 



where A345 = A3 -|- A4 -|- A5 . 

When testing a model such as the CMSSM, one imposes two classes of constraints: theo- 
retical and experimental. Let us discuss them. Minimizing this scalar potential leads to two 
theoretical constraints given below and their solutions for 7;,tan^ are the same as those of 
the eqs /i = /2 = where /i and /2 are introduced for later convenience: 



,2 , ^2, 



1/2 



5/3 """ 5:5 



C/TY) (7 \ 

2A^— -777^-- =0, /2(7i;7;,/3,yt,y6,---) = tan/3-tan/3o(7i,7;,yi,y6...), (5) 



Here v = \Jv\ + v\ is a combination of vev's of /7]', l^- The order of the arguments of /1.2 is 
relevant later, while the dots denote other parameters (other Yukawa or gauge couplings,...) 
present at one-loop and beyond, that we ignore in this section only, without loss of generality. 
As a result of these two constraints, the EW minimum solutions for v and tan/3 become 
functions of the (mass dimension 1) parameters 7^ of the model which for CMSSM are shown 



above, in a standard notation. When discussing fine tuning, usually only the first constraint 
in (0) is considered, although the second is equally important, as our result for p{D) will 
show. These constraints fixing v, tan/3 are assumed to be factorized out from the likelihood 
function C{D\ji) and can be imposed by Dirac delta functions of arguments: 

Hfi{7uV,P,yt,yb)), S{f2{-fuV,/3,yt,yb)), i = T75. (6) 

There is also a second class of constraints, that comes from the experiment, such as the 
accurate measurement of the masses of the top {mt), bottom (mh) and Z boson (mz)- Given 
the high accuracy of these measurements, one can assume some Gaussian distributions for 
the associated priors when evaluating the probability density T'{'yi\D) or the evidence p{D). 
However, for a more qualitative analysis and to good approximation one can again implement 
these constraints (likelihood) by Dirac delta functions of suitable arguments 

5{mt-m'^); 5{mb - m'^); 5{mz-m%), (7) 

where m^, m^,m^ are experimental values. One can consider similar constraints for aem and 
as gauge couplings but for simplicity we do not do that (their implementation is similar). 

When testing the SUSY models with a given set of parameters (such as 71,. ..5 for CMSSM), 
one should in principle marginalize (i.e. integrate) the density ^(74 1-0) over the "nuisance" 
parameters. Examples of these nuisance parameters are those already present in the Standard 
Model. These are the Yukawa couplings yt,yb, ■■■ OE] which were restricted (in the analysis of 
this section only) to the simpler case of top and bottom Yukawa couplings. Other parameters 
to integrate over are the dependent parameters: the vev v and tan/3 which can (in principle) 
take any value, until fixed by minimization constraints ([5]), ([6]), also ([7]) for v. 

To compare two SUSY models, one should compare their evidence p{D) for similar priors 
and data D. To compute p{D), one integrates over all parameters (of the SM and those 
mentioned above) and over 7, as well, with chosen priors p(7i). For the CMSSM case, after 
imposing the above constraints with the corresponding priors, one finds 



p{D) = N dji...dj5 p(7i, ...75) dyt dy^ dv (i(tan/3) p{yt) p{yb) 

X 6{mz — iT^z) ^('^t — "rrit) ^{"n^b ~ ™^b) 

X 5{fi{li\v,l3,yt,yb)) 5{f2{ii;v,l5,yt,yb)) C{D\-fi,2,...5;P,v,yt,yb)- (8) 

where C{D\ji; 13, v, yt,yb--) is the likelihood of fitting the given data (D) with a particular set 
of values 7^; i=l,...5, etc; the priors p(7i,...75) and p{yt,b) are not known, but logarithmic 
or flat priors are common choices for individual parameters. Regarding the priors p{v) and 
p(tan/3), these are already included and given by the corresponding Dirac delta's shown in 
(©5 ©• We integrated over yt,yb rather than over the corresponding masses mt,mb. This is 
a possible choice, preferable because the masses are derived quantities, see discussion in [3]. 
Finally, leaving aside the integral over jj and p(7i, •••75), the above equation simply gives the 
probability density 1^(^110). 



The important point about eq.® is that now ah parameters 7,, v, tan /3, yt, 2/b, • " " that we 
integrated over can be regarded as arbitrary, since the constraints that render them dependent 
variables are implemented by the Dirac delta functions associated to the theoretical and 
experimental constraints. £ is a function of the CMSSM parameters, but also of the nuisance 
parameters {yt.b) and u,tan/3. Finally A^ is a normalization constant not important below. 

To evaluate p{D) , one uses mz = g v/2, rrit = ytv sin /5/\/2, mi, = ytv cos /3/\/2 and after 
performing the integrals over yt, yf, and v one finds 

p{D) = -^ / c?7i....c?75 p(7i, ...75) d(tan/3) p{yt{P)) p{yb{P)) csc(2/3) 
X 5[/i(7i; /3, ^0, yt(/3), ym)] S[f2{7^; (3, vq, yt{P), Ul^))] 
X £(D|7i;/3, 7;o,yt(/3),yfe(/3)), (9) 

with g"^ = Qi + 52 where gi (52) is the gauge couplings of U(l) (SU(2)) and 

fo = 2?4/5 = 246GeV, UP) = V2 m^J (vq sin /3) , UP) = ^m^b/i^ocos 13). (10) 

Integrating ovei|f| /3: 

piD) = ^ f dji....droPili,-To) {p{UP))p{UP))) csci2l3) [(12)'^]-' 



X 4/1(7,; /3, vo, UP), UP))] ^{D\^uP, vo, UP), VbiP))}^ 
dS^, 71—75 pill, ■■■l5)[p{UP)) p{yb{P))) csc(2/3) 



9V0 Jm 

X [if2)'p\V^^nviJf,(3oil^))\]~' C{D\^,;P, vo, UP), UP))} ^ ^, ■ (H) 

Above (/2)« denotes the partial derivative wrt the variable tan/3 of the function /2 of argu- 
ments: f2{li; P,vo,yt{P),yb{P)), where ([TO]) is used. The curly bracket is evaluated at the 
unique root /3 = /3o(7i) of the second minimum condition in ([5]) of the scalar potential: /2 = 0. 
Through this condition, /3 becomes a function of the independent parameters 7,, as usual (one 
can eventually trade /3q for Bq, as often done, but we do not do this here). In the last step we 
converted the integral into a surface integral where Ai is the surface defined by the equation 
/i = while dS-y is the surface element in the parameter space {ln7j}. Recall that /i = is 
one minimum condition which together with the second one f3 = Poili) (or /2 = 0) control 
the value of p(Z?). A notation was used V^Jiin; P,vo,UP),yh{P)) = V^.w(7i;/3o(7i)) where 
V = —m?/X has the arguments shown and V^. is the gradient in coordinate space {ln7j}. 



■^We use S{g(x)) = S(x ~ xo)/\g'\ _ with g' the derivative wrt x evaluated in xq; xq is the unique root of 
g{xo) = 0; we apply this to a function g{j3) = 72(7;; /3, vq, yt{l3),yb{l3))) for x = tan/3 with the root /3o = Poi'Ji)- 

"One uses f^„ f{zi, ..., Zn) S{g{zi, ..., Zn)) dzi....dzn= fg^_^ /(zi, ...Zn) dSn-il/IV^.pl with Zi -^ Inji where 
Sn-i is defined by <? = and V is in basis Zi = ln7i. Another form of Hll|) is found by replacing dS, V by 
their values in {7^} space (instead of {ln7i}) and removing the product 71. ...75 in integral (|ll|l . 



The important result is that p{D) contains a suppression factor 1/A.q where we denoted 

=> p{D) r^ dS^ — C X {priors), (12) 



with 

A fS^A^V^'^ A dlnv{"fk;/3o{lk)) _ . rj /1Q^ 

\= [Z^ \j) ' A^^. = ^^-— ; -fj = mo,mi/2,fio,Ao,Bo. (13) 

Note that Aq(7j) contains a derivative of v ~ /i evaluated at /? = /3o(7j), so it encodes 
the effects of variations about the ground state of both minimum conditions ([5]), see the two 
Dirac 5's in Q. A good stability of these conditions under such (quantum) variations requires 
small Aq. Interestingly we also notice that Aq(7j) = Ag so A^ is preferred by the calculation 
of the Bayesian evidence p{D). The points {7j} of smaller A^, give larger contribution to 
p{D), but this also depends on C or priors. We can say that 1/Ag is an extra effective prior, 
emerging when marginalizing over parameters, subject to the theoretical constraints. With 
p{D) ~ 1/Ag, points of large Ag pay the fine-tuning cost and so have a small impact on p{D). 
The latter is then used to decide which of two models is more probable. 

For further illustration, assume log priors for Yukawa couplings p{yt{b)) = ^/Vtib) ^'^^ foi' 
SUSY parameters {7j}, using p(7i...75) = p(7i)---p(75) and with ^(7,) = 1/7j. Then 

To conclude, 1/Ag is an extra effective prior Peffiji) of the model and Ag emerges as a 
measure of fine tuning. In the general case Peffi'fi) can be read from (jlip . ()12p and the link 
between Ag and the Bayesian evidence p{D) is clear. Numerical studies of p{D) or £(Z)|7j) 
should then include such effect due to the two theoretical constraints. To our knowledge 
this effect was so far overlooked in such studies. 

Note that Ag that emerges in eqs. (jl2p . (|14p . does not contain partial derivatives wrt 
Yukawa couplings. This is because these are nuisance (SM-like) parameters that were inte- 
grated out, so are included as a global effect. Also, such parameters are not part of the new 
ones (7j) that SUSY introduces, so it is no surprise that they are not explicitly manifest in 
p{D) or in the denominator under integrals (|12p . (|14p . 

These above results bring technical support to a physical meaning of the fine tuning. They 
show that it is desirable to have a smaller A, as also expected from physical considerations. 
Again, one should remember that this may not always be the region from where p{D) receives 
the largest contributions, as this depends also on the priors, the integral(s) or their measure. 
Note also that changing the priors of the nuisance parameters or the measure can give different 
values for Bayesian p{D) although with enough data D one expects this dependence to become 
weaker. With this technical motivation for the fine tuning measures and their relation to p{D) , 
C{D\ji) below we study the values of Ag, A^ax for many SUSY models. 
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Figures 1 to 4: Ag, A^ax versus M^iggs', lightest grey (0) area: excluded by SUSY mass bounds; 
darker grey (1): excluded by 6— >-S7, B ^ fi^fj,~, Sp; dark grey (2): excluded by condition 5a^ > 0. 
Coloured area: allowed by data and Sufj, < (25.5 + 2x8)10"^"; 5a™°^ is shown colour encoded. Area 
outside contour: Sa"^'^^ < (25.5 — 2x8)10~^° (2ct). Red area (inside): largest Sa/^ is within 2ct of 6a'^^^. 



3 Numerical results for A in generic supersymmetric models. 

We present our numerical results for Ag and A^ax in a comparative analysis for generic 
models used for SUSY searches at the LHC. We scan the entire parameter space {7j} of the 
models, consistent with the theoretical constraints, using a two- loop leading- log analysis. Ag 
and Amax are presented as functions of physical scales (mass of higgs, stop, gluino, SUSY 
scale msusy = {iTT-i "^i,)^ ) with constraints (muon magnetic momentum da^, dark matter 
abundance ^h?). The models considered are: 

• the CMSSM model, of parameters 7j = {mo,m-i/2,^0) ^O) ^o}- Then Aq is that shown in 
([!]) and (fTBj) . evaluated at the two-loop leading log level. See [3] for a recent study, whose 
results were recovered by this work. 

• the NUHMl model: this is a CMSSM-like model but with Higgs masses in the ultraviolet 
(uv) different from ttiq, m^^ = m^^ ^ niQ, with parameters 7j = {mQ,mi/2, IJ-Oj^o, BQ,m^^}. 
Then Aq is as in (J13p with summation over this set. 

• the NUHM2 model: this is a CMSSM-like model with non-universal Higgs mass, rii™ ^ 
m'^^j^mQ, with independent parameters jj = {rno,ini/2, (J-OjAq, BQ,7n^^,m^^J^}. Then Aq is 
that of (jl3|) with summation over this set. 

• the NUGM model: this is a CMSSM-like model with non-universal gaugino masses m\., 
i = 1,2,3, with 7j = {niQ, fiQ,AQ, BQ,mx^, mx^^mx.^}. Then Aq is given by ([T3|) with the 
sum over this set. 

• the NUGMd model: this is a special case of NUGM-like model with a relation among the 
gaugino masses mx^, i = 1,2,3, of the type mx^ = f/iWii/2) where T/1,2,3 take only discrete, 
fixed values. Such relations can exist due to some GUT symmetries, like SU(5), SO(IO), 
etc. The particular relation we consider is a benchmark point of [8] with mx^ = {l/3)mi/2, 
mxi = (—5/3)7711/2) f^\2 = "^1/2) corresponding to a particular GUT (SU(5)) model, see 
Table 2 in [8]. As a result, Aq is that of ([T3|) with 7^ = {mo,mii2-, Aq, Bq, ^q} . 

In all models we also evaluate the alternative definition of A given by 

Amax = max I A^ 1 , 7: parameters of mass dimension 1. (15) 

and where the set 7^ is listed above for each model. 

Before presenting our results let us describe the method used. The scan over the full 
phase space of each model was done using Pythia 8 |14] random number generator. The 
public code micrOMEGAs 2.4.5 [15] and SoftSusy 3.2.4 code [16] were then used, with the 
latter adapted to compute for all models the fine tuning of the electroweak scale at the two- 
loop leading log level (instead of its default, one-loop calculation). This includes two- loop 
tadpoles to the two electroweak minimum conditions. The data output was then filtered by 
the experimental constraints. The run time to generate the phase space points of the five 
models was about 15000 one-day jobs on the CERN computing service, and each plot was 
generated from ~ 4 x 10^ points in a random scan of the parameter space (for alternative 
and recent data analysis see [T7]). 
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Figures 5 to 8: A,, Amax versus Mhiggs- Various grey areas and Sa^ values: as for Figures 1 to 4. 
Colored ares: allowed by data other than Sa^. Blue area: flh^ < 0.1099 — 3 x 0.0062. Red area 
0.1099-3 X 0.0062 < H/i^ < 0.1099 + 3 x 0.0062 (3ct saturation). Yellow: nh^ > 0.1099 + 3 x 0.0062. 



Experimental constraints 


Values used 


SUSY particle masses 


Routine in micrOmegas 2.4.5, "MSSM/masslim.c" 


Muon magnetic moment 


6a^ = (25.5 ± 2 X 8) X 10~^° at 2o- [IB]. 


6—7-8 7 process 


3.03 < 10^ Br (6 -^ S7) < 4.07 at 2o- p]. 


Bg —7- ji'^ ^^ process 


Br(5, -^ fi+fi-) < 1.08 X 10-s at 2a [20]. 


p-parameter 


-0.0007 < Jp < 0.0033 at 2a [21]. 


Dark matter relic density 


rj/i2 = 0.1099 ± 3 X 0.0062 at 3cj [22]. 



Table 1: Experimental data constraints. (5a^ includes the theoretical error and is not imposed on the 
data, but its values are shown as a contour plot (at 2a) or colour encoded from which larger deviations 
can be read (Sa). For the other processes in the table, only the experimental error is considered, and 
the details of their theoretical calculation are provided by micrOMEGAs 2.4.5 [T5|, see also its manual 



for v.2.4 available at http://lapth.in2p3.fr/micromegas/. The central values for nitop = 173.1 GeV 



and a3(m2) = 0.1184 [21] were used as inputs in SOFTSUSY. Note that a combined la increase of 
top mass and la decrease of a3{m%) can decrease A^ax by a factor as large as 2 (best case scenario), 
see later. 

Our resultcl shown in the plots allow the reader to set his own constraints on physical 
scales such as the liiggs mass, gluino, stop mass or SUSY scale rugusy, ^cl^ or dark matter 
abundance and infer from that the amount of fine tuning. Note also that the LEP2 bound on 
Mhiggs is never imposed on our figures, and we let the reader to do this, in the light of future 
LHC resultqj. This has the great advantage that the impact of future bounds from LHC 
on these physical scales can very easily be seen on the plots, without the need to re-do the 
whole analysis. Space constraints do not allow us to also present a description of the allowed 
parameter space {'ji} = {Aq, BQ,etc..} used in these plots, due to complicated correlations 
among these, that can only be presented as more additional figures, that we postpone to a 
future work. Finally, the parameter space ({7j}) that we scanned over was: ^o £ [—7, 7] TeV, 
mo e [0.05,5] TeV, mi/2 G [0.05,5] TeV and also 2 < tan/3 < 62. All plots are marginalized 
over tan/3 and {ji}. 

3.1 A versus Mhiggs and the values of 5a^. 

In figs. 1 to 4 a), b) and 17 a), b), we show the plots for A^ and A^ax as functions of the 
mass of the lightest Higgs boson Mhiggs, for all models: CMSSM, NUHMl, NUHM2, NUGM, 
NUGMd. The impact of (5a^ constraint is also shown with a contour line displaying an island 
of its largest values, within 2a of 6a^^. For other values {3a deviations, etc), the largest 5a^ is 
also shown colour encoded, see the scale on the right side of the plots. The lightest grey (level 
0) areas in these plots are excluded by the lower bounds on the spartners masses obtained 
from negative SUSY searches. The darker grey (level 1) areas are excluded by Bg^ fi~^fi~ , 



^ After this work was completed, an updated bound on Bs — > /i'^yU was published [23]. We checked that 

our fine tuning estimates are unchanged, for a higgs mass in the now preferred region of 122 GeV to 128 GeV. 

^ Note that a flat bound like LEP2 bound on Mhiggs should be used with care since it applies only to SM. 
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Figures 9 to 12: A, versus gluino (left) and SUSY scale (right) for various models; Mhiggs < 123 
GeV in area below the dotted line; similar plots exist for Amax', various grey areas as in Figures 
1 to 4, forbidden by data. Black area: Mhiggs < 111-4 or Mhiggs > 130 GeV. Outside the red area: 
g^max < (25 - 2 X 8)10^1°; inside this area Sa]^"-" is within 2a of Sa^^^'P. See also caption of Figs 1-4. 
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b^ sj and 6p constraints. The dark grey area (level 2) that we also show, visible only for 
NUGM model corresponds to (5a^ < and has (mx^fi) < and is present at M^iggs < 115 
GeV. This region is excluded by demanding 5a^>0, ((ttiaj/u) >0), preferred by (5a^^^ data. 

As it is shown in these figures, the LEP2 bound (114.4 GeV) |24j on the higgs mass was 
not imposed. Note however that above this value both A's are largely independent of the 
experimental data (ignoring (5a^) for all models other than NUHMl, NUHM2; these still have 
some dependence on data (the small grey area). This is interesting and suggests that the 
range of values of Ag, A^ax can be fixed mainly by theory and the higgs mass bound, with 
little or no impact from other data. 

As it is seen from these results, the differences between Ag and A^ax are practically 
negligible. For a given model and a fixed Higgs mass, there is a relative factor between 1 and 
2 and which can be safely ignoreqI|. There is also very similar behavior i.e. various contour 
lines such as that of maximal 5a ^ are nearly identical for both Ag, Amax- This is interesting 
and shows that one can use either definition for fine tuning to obtain a rather similar result. 

In the CMSSM the minimal value of Ag and Amax is situated for higgs mass near the 
LEP2 value as also discussed in |5]. This means that to respect the LEP2 bound on the 
Higgs mass there is no fine tuning cost due to quantum corrections. This corrects common 
but wrong opposite claims in the literature. Further, if one accepted the principle that A of 
a model should actually be minimized, then one immediately has a CMSSM prediction for 
Mhiggss^ 115 ±3 GeV without using experimental constraints (ignoring here 5a^), for details 
see [SJf|. In models other than CMSSM and after imposing the LEP2 bound, the fine tuning is 
again smallest near this scale. For Mhiggs near 115 GeV, Ag ^ A^ax ~ 10 to 100, depending 
on the model. Above this mass value, both Ag, A^ax grow very fast (~ exponentially), 
due to the quantum corrections to the Higgs mass (which are logarithmic in mgusy)- As a 
result, for the currently interesting region discussed by CMS and Atlas experiments [9l|10], of 
123 < TTi/j < 127 GeV, there is significant amount of fine tuning required, Ag ~ A^ax ~ 200 
to 2000; for Mhiggs = 125 GeV, Ag w A^ax ~ 500 to 1000, depending on the model. From 
these results one could say that NUGM is preferable also because it could more easily comply 
with da^f^^ (2o")- Finally, let us mention that a combined la increase oinitop and la reduction 
of a3(m^) can reduce (best case scenario) the fine tuning for a fixed higgs mass by a factor 
near ~ 2 or so for the CMSSM [5], with similar effect expected for other models. 

3.2 A versus Mhiggs and dark matter relic density. 

Let us now discuss the results of figures 5-8 a), b) and 18 a), b). These present the impact 
of the dark matter relic density constraint. Again, no significant difference between Ag and 
Amax is observed for the models considered. The meaning of light and dark grey areas is 
the same as in the previous figures. In blue we show points that are consistent with dark 
matter relic density within 3a, i.e. these points have Qh^ < 0.1099 — 3 x 0.0062. The red 
points saturate the relic density within 3a deviation from the central value. Finally, yellow 



^ In general no individual A^^ dominates clearly for all higgs masses, see fig.2 in [5] for the CMSSM. 
*There is a ±(2 to 3) GeV theoretical uncertainty from the various public codes [TS] [551 [5S] . 
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Figures 13 to 16: A^ versus Mhiggs with largest (left plots) and lowest (right plots) gluino mass; For 
given Aq, M^iggs one infers the gluino mass range. Area allowed by data (except (5a^) as shown above 
the continuous line, see also Figures 1-4. Values as large as Sa^ < (25 — 2 x 8)10~^° are outside the 
closed contour; inside the contour: largest Sa^ is within 2a of Sa'^^P and the gluino mass satisfies it. 
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points correspond to a relic density larger than that of the red points. Notice that with the 
exception of the CMSSM case at M^iggs ~ 115 GeV region, for a higgs mass above this value 
one can easily saturate the relic density. This is true in particular for points near the 125 
GeV region, although in CMSSM this may become more problematic (too large Qh'^). 

The contour area of maximal values of 5a^ already shown in previous figures 1-4, 17 is 
also presented. However, it should be stressed that points inside this contour area that have 
the dark matter abundance as shown (in blue, red or yellow) are not necessarily the same 
points that also have the largest 6a ^ within 2a of 6a^^\ (these are different projections on 
the 2D plane shown) . However points that satisfy a relic density constraint and also have the 
largest (5a^ within 2cr of experimental value do fall within a smaller area inside the contour 
shown. Of all models, NUGM could both saturate the relic density and fall within the 5a^ 
contour line, for large range of higgs mass, although the fine tuning cost grows exponentially 
with Mhiggs- Figures 5a), b) and 8 a), b) show again that CMSSM and NUGM are the least 
sensitive models to any experimental constraints other than 5a^, for Mhiggs larger than « 115 
GeV (negligible grey areas). Finally, since Ag, Amax are so similar, below we shall present 
only results for Ag. 

3.3 A versus Mgiumo and A versus nisusy 

So far we investigated the fine tuning as a function of the higgs mass. However, it is useful 
to present its dependence on other particles masses, and we do this for the gluino and the 
SUSY scale nisusy. This is useful since LHC searches for gluino or other SUSY partners can 
have a strong impact on fine tuning. This is seen in figures 9-12 a), b), and figures 19 a), b), 
where we show the dependence of Ag on the gluino mass (figures a)) and on rususy (figures b)) 
for all models. The light and dark grey areas have the same meaning as before, while the 
areas in black are ruled out by the higgs mass constraint 111.4 < Mhiggs < 130 GeV that we 
imposed (this allows 2-3 GeV uncertainty for Mhiggs at two- loop leading log level [TBI [25l [26] ) . 
Contour (dotted) lines of a maximal value of 123 GeV of Mhiggs are displayed for all models: 
the points below this line respect this bound while those above can have larger values. The 
advantage of these plots is that if future data rules out Mhiggs < 123 GeV, the whole region 
below (outside) the dotted line (contour) will be removed from the plots, to leave a small, 
restrictive region. 

5a^ is also shown in colour encoded areas, with a red island area showing the largest 
possible value with 5a'^"'^ within 2a of the experimental central value. Note again that the 
(5a^ contour and the dotted line of upper bounds on higgs mass are different projections on 
the 2D plane of the figures. That means that points that have largest 5a^ within 2a of 5a^f^^ 
are not necessarily the same points that simultaneously have Mhiggs as large as 123 GeV. 
The impact of future gluino mass or nisusy bounds from the LHC can easily be seen on these 
plots, together with the associated fine tuning cost. The models NUGM and NUGMd relax 
the lowest bound on the gluino mass due to their non-universal gaugino masses. 
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Figures 17-20: The benchmark NUGMd modch the description of the plots is identical to that in Figs. 
1-17; but applied to NUGMd model, as follows: Figs 17a), b) - as for figure la),b). Figs 18a), b) - as 
for figure 2a),b). Figs 19 a), b) - as for figure 3a),b). Figs 20a), b) - as for figure 4a),b). 
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3.4 A versus Mhiggs and the gluino mass range. 

A complementary presentation of the results of figures 1-4 a), 17 a) and figures 9-12 a), 19 a), 
b) is that of figures 13-16 a), b) and 20 a), b). In these Ag is presented again as a function 
of the higgs mass, but with the gluino mass as a parameter, with its largest value in plots 
a) and lowest possible value in plots b), see the colour encoded scale. In this way one has 
a clear picture of the whole range of allowed values of gluino mass for a given A and higgs 
mass. Intermediate values of gluino mass are colour encoded. For the large Mhiggs, above 
125 GeV the gluino mass tends to be larger (above 1 TeV), and within a narrow range, with 
increasing fine tuning cost. The range of values of gluino is rather similar in CMSSM, NUGM 
or in NUHMl, NUHM2. 

One important remark about the contour of largest 5a^ shown: the gluino mass range 
shown inside this contour respects all experimental constraints, including the constraint of 
(5a^ (within 2a)\ If this constraint is lifted, the range of gluino mass, for a fixed higgs mass 
and fine tuning, would be larger. This also explains the sudden change of colour /spectra of 
gluino masses around the contour line of 5a^ as compared to region immediately outside the 
contour. 

3.5 Stop versus gluino mass, with the largest Mhiggs and minimal A. 

For a future comparison with results from LHC searches for new physics, we also present in 
figures 21, the dependence stop versus gluino mass and with the largest value of Higgs mass 
that is possible with the former two fixed. The minimal fine tuning cost that comes with this 
is also shown in the corresponding areas (bordered by red contour lines), while the largest 
Mhiggs allowed is colour encoded, see the scale on the right side of the plots. The latest 
bounds on the gluino and stop masses can be translated into (upper) bounds for the higgs 
mass. Currently, stop-gluino exclusion plots from the LHC exist only for simple models that 
cannot be used for comparison ^7\, see the first plot in figure 21. Eventually, at very large 
gluino and stop masses the (minimal) fine tuning cost becomes too large and the models may 
be considered unrealistic. It can be seen from these plots how the lowest allowed fine tuning 
increases as the higgs mass goes towards its upper limit. If one rules out values of fine tuning 
of say Aq > 100 one immediately removes the area outside the contour line that borders this 
region, to leave a significantly smaller area of correlation stop-gluino-higgs mass. 

While the CMSSM, NUHMl, NUHM2 are more restricted by superpartners masses (ex- 
cluded light grey areas), in the case of the NUGM and NUGMd, not surprisingly, the impact 
of the spartners mass bounds is small (since the universality condition was relaxed). Again, 
the NUGM model is less restricted, allowing a large higgs mass (125-128 GeV), with a stop 
as light as 400 — 500 GeV and gluino mass between 2 — 3 TeV. 

4 Final remarks and conclusions 

Low energy (TeV scale) supersymmetry is thought to solve the hierarchy problem without 
undue amount of fine tuning (A). However there are different opinions on what the best 
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Figure 21: Top left plot: Atlas stop-gluino exclusion limits in a simple supersymmetric model |27| . 
Rest of the plots: The dependence stop vs gluino mass in CMSSM, NUHMl, NUHM2, NUGM, 
NUGMd models, in this order. We present the lowest value of Ag in the areas where it is shown, 
bounded by red contour lines and with no upper bound. Areas of largest Higgs mass are also shown, 
colour encoded, see the scale on the right side (minimal value: 111.4 GeV). One can easily see the 
largest higgs mass and the minimal fine tuning cost, for given gluino and stop masses. Grey area is 
excluded by SUSY mass bounds. Black area is excluded by imposing the constraint Mhtggs > 111.4 
GeV and M^iggs < 130 GeV. This dependence can eventually be compared with similar future plots 
from CMS/ATLAS searches. 



definition for A is, or what upper value is allowed for it while still claiming a SUSY solution 
to this problem. To avoid a subjective choice on these two issues, we performed a study of 
A using two common definitions Amax and Ag and made no assumption about their largest 
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allowed values. We also discussed the relation of A to global probabilities (in the parameter 
space) to fit the data. We analyzed generic models: CMSSM, NUHMl, NUHM2, NUGM 
and a benchmark model, NUGMd, at two-loop leading log level, and both A^ax and Ag were 
presented as functions of the higgs, gluino, stop mass or the SUSY scale, with additional 
constraints like dark matter or 5a^ . The advantage of this comparative analysis is that using 
the figures for Ag, Amax, future experimental constraints can immediately be converted into 
an updated estimate for the fine tuning level of these models, without the need to re-do the 
whole analysis. The reader will then decide whether the amount of tuning so obtained is still 
acceptable for a solution to the hierarchy problem. 

The measures of fine tuning were originally introduced more on physical intuition than 
rigorous mathematical grounds. In this work we provided mathematical support for the 
fine tuning via a quantitative relation to Bayesian evidence p{D). As direct result of two 
theoretical constraints (EW min conditions) , we showed that a fine tuning measure Ag = Ag 
emerges as an additional suppression factor (effective prior) of the averaged likelihood under 
the initial priors, under the integral of global probability of measuring the data (the evidence 
p{D)). So the Bayesian evidence calculation prefers Ag as a fine tuning measure. As a result, 
the evidence p{D) ~ 1/Ag therefore points of large Ag (strongly fine tuned) have little or no 
impact on the global probability of the model to fit the data. These results provide technical 
support to the idea that fine tuning has a physical meaning and that preferably it should 
have small values in realistic models for the corresponding point in the parameter space. 

Our numerical results for Ag and Amax as functions of the higgs mass, showed that they 
have close values for the same higgs mass and also very similar behaviour for all models 
considered. There is a small discrepancy factor between them (between 1 and 2) which 
is most visible for regions of the higgs mass that are anyway excluded by the data. All 
these results show a good independence on the actual definition used for fine tuning. For 
115 < Mfiiggg < 128 GeV there is a relative independence of Ag or Amax on the experimental 
constraints (other than 5a^) for CMSSM or NUGM, NUGMd, with a minor dependence for 
NUHMl, NUHM2. So in this case Ag, A^ax are largely controlled by theoretical constraints. 
Also, the dark matter relic density can in all cases be saturated within 3a of the current 
value. 

The dependence of both A on the gluino mass or on the SUSY scale shows a similar 
behaviour for all models. The CMSSM, NUGM and NUGMd models show a lower amount 
of fine tuning for the same experimental constraints, and NUGM can even accommodate 5a^ 
(2(t) and Mhiggs ~ 125 GeV, however in this case there is always a fine tuning cost. As our 
plots showed, for the CMSSM no fine tuning amount can reconcile 5a^ (2(t) values considered 
(i.e. contour at 2a in the plots), with a Mhiggs > 120 GeV region which is situated outside 
this contour. For a Higgs mass near 125 GeV, the fine tuning is of order 0(1000) in all 
models other than NUGM, NUGMd where it is of order 0(500). There is a strong (roughly 
exponential) variation of A with Mhiggs . A reduction of 2 GeV of Mhiggs can bring down both 
A's to A ?» 200 to 500, depending on the model. For Mhiggs = 115 GeV, Ag ~ Amax ~ 10 
to 100 and in the CMSSM this A corresponds to a global minimum. Finally, let us mention 
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that the combined effect of a la increase of the top mass and a similar reduction of the 
measured strong couphng at EW scale can reduce the fine tuning for a given M^iggg by a 
factor near 2 or so in the CMSSM case [5]. Although we did not studied it here (due to long 
CPU runs), we expect similar effect for the other models. This is because Yukawa corrections 
help radiative EW breaking (reducing A) while QCD corrections have the opposite effect in 
the loop diagrams. 

Are the values of fine tuning that we found too large? Based on previously agreed but 
highly subjective "reasonable" values of A ~ 10 — 100, the answer is probably affirmative. 
However, a clear answer is difficult, largely because A depends ~exponentially on the higgs 
mass, so any small correction to it has a strong impact on A. But comparing all models, 
for the same experimental constraints, there seems to be a preference for NUGM case when 
also considering the (5a^ constraint. We let the reader to make his own opinion, based on 
the above results and figures and also on future LHC data (on gluino, higgs, stop and nisusy) 
whose updated impact on our A can easily be obtained. Also it should be kept in mind that 
very simple new physics beyond these SUSY models (like CMSSM with a gauge singlet with a 
TeV-scale SUSY mass term or a massive U{1)') can lead to a very acceptable A ~ 0(10) for a 
higgs mass as large as 130 GeV |llj . Further, subjective criteria also exist in other approaches 
that compare the probability of various models, such as those based on the Bayesian approach. 
Indeed, the evidence p{D) also has some dependence on the priors choice (flat, log, etc), until 
eventually more data can improve our knowledge of the models. We hope that the clear link 
between fine tuning A^ and p{D) that we established together with our plots for both A's 
will provide the starting point of a more detailed study. 
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